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Abstract 



We performed the observation of the flux densities of Sgr A* at 90 and 102 GHz in order to detect the 
time lag between these frequencies using the Nobcyama Millimeter Array, which was previously reported at 
lower frequencies. We detected a radio flare during the observation period on 6 April 2005 and calculated the 
z-transformed discrete correlation function between the light curves. The time lag between these frequencies 
was not detected. If the expanding plasma model which explains the time lag at lower frequencies is valid, 
t— » \ the light curve at 90 GHz would be delayed with respect to the one at 102 GHz. This result suggests that 

the plasma blobs ejected near the Galactic Center black hole may be widely diverse especially in optical 
thickness. Another possibility is that the major portion of the flux above 100 GHz docs not originate from 
the blobs. 

Key words: black hole physics — Galaxy: center — galaxies: active — galaxies: nuclei 



Introduction plasma ejected near GCBH. The Very Large Array (VLA) 

detected time lags of 20-30 min between flare peaks at 22 
Sagittarius A* (Sgr A*) is a compact source with emis- and 43 GHz (Yusef-Zadeh et al. 2006; Yusef-Zadeh et al. 



sions from radio to X-ray, and it is believed to be as- 2008). These time lags strongly suggest that the radio flux 



Z^ [ sociated with the Galactic Center black hole (GCBH). is emitted from adiabatically expanding plasma. Several 

Sgr A* is considered a laboratory for studying the activ- authors have successfully applied the model of adiabat- 

ity of galactic nuclei. In the radio regime, very long base- ically expanding plasma to interpret the light curves of 

line interferometer (VLBI) is a powerful tool for exploring the observed IDVs (Yusef-Zadeh et al. 2008; Eckart et al. 

the structures of galactic nuclei. However, VLBI obser- 2008; Eckart et al. 2009; Li et al. 2009). However, for the 

vations of Sgr A*, at least up to 40 GHz, are not fruitful determination of the time lag, truly simultaneous obser- 

in spite of great efforts because the intrinsic structure of vations at both frequencies were performed in only one 

Sgr A* is hidden by electron scattering, as indicated by epoch. For verifying such observations, additional simul- 

the observed size of Sgr A* obeying the inverse square taneous observations are needed, preferably at several dif- 

law at those frequencies (e.g., Shcn et al. 2005). The ferent frequencies to provide further constraints on source 

VLBI technique at submillimctcr wavelengths is advanc- structure of Sgr A*. However, if the separation in obscr- 

ing rapidly (Doeleman et al. 2008), but has yet to obtain vation frequencies is too wide, it may lead to ambiguity 

an image of Sgr A*. On the other hand, flux monitoring in identifying true counterparts among IDV at each fre- 

of Sgr A* at short centimeter and millimeter wavelengths quency. Thus, we made observations at 90 and 102 GHz 

has been useful to explore the region around GCBH. Intra- using the Nobeyama Millimeter Array (NMA), Nobeyama 

Day Variability (IDV), flux variability with a timescale of Radio Observatory (NRO) 1 to measure the time lag in ra- 

a few hours, has been observed frequently (e.g., Miyazaki dio variability between these frequencies. NMA can ob- 

ct al. 2004; Mauerhan et al. 2005; Li et al. 2009). The serve Sgr A* at these two frequencies simultaneously, and 

emitter of the radio flux of Sgr A* is in controversy but un- thus should be able to robustly detect the time lag even 

doubtedly the emissions come from very near the GCBH. 

rr, .,.,.,., , 1,1 t .j j NRO is a branch of the National Astronomical Observatory of 

1 wo possibilities nave been proposed as the radio emitter: , m- *• i t an. n, <_ ,o ■ 

*l i Japan, Nationai fnstitutcs of JNaturai sciences. 

hot plasma in the accretion disk of GCBH or expanding 
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if it is relatively small. 

2. Observation and Light Curves 

Aiming to detect a small time lag, we used NMA on 6 
April 2005 under stable clear weather to observe the flux 
densities of Sgr A* at 90 and 102 GHz. NMA consists of 
six 10 m element antennas equipped with double side band 
(DSB) superconductor-insulator-superconductor (SIS) re- 
ceivers with a single linear polarization feed. The Ultra- 
Wide-Band Correlator (UWBC), which has 1 GHz band- 
width, was employed for the backend (Okumura ct al. 
2000). The lower and upper side band signals (90 ±0.5 
and 102 ± 0.5 GHz) were separated by 90° phase switch- 
ing to obtain simultaneous data. The instrumental gain 
and phase stability were calibrated by alternating obser- 
vation of Sgr A* and NRAO530 (2.16 Jy in April 2005) 
at scan intervals of 23 min. The observation of Sgr A* 
was performed in the NMA array configuration with inter- 
mediate baselines, "C-configuration" , giving a projected 
baseline range of ~ 4 — 47 fcA. The UV data were pro- 
cessed with the UVPROC-II software package developed 
at NRO (Tsutsumi et al. 1997). To obtain the flux den- 
sities averaged every 1 min for a light curve, we applied 
a point source model to the visibility data, which was 
self-calibrated for phase, by using the MIRIAD package. 
The visibility data are restricted to the projected baseline 
larger than 25 fcA in order to suppress the contamination 
from the extended components around Sgr A*. The ab- 
solute uncertainty of the flux scaling is about 10% in the 
100 GHz band. In addition, the phase noise due to at- 
mospheric fluctuations, which dominantly contributes to 
the visibility fitting error, and the effect of contamination 
from the extended components are considered. These er- 
rors, which are estimated to be order of a few percent, 
are smaller than the flux scaling error. Thus the relative 
uncertainty of the flux density in adjacent bins should be 
lower than the absolute uncertainty. The signals at 90 and 
102 GHz, which are simultaneously received by the DSB 
receiver, arc separated by 90° phase switching in the local 
oscillator as mentioned above. Most of signal paths at the 
both bands are the same, for example the same antenna 
and the same receiver system, and thus the gain drifts of 
the both bands are almost common. At the NMA site, 
the atmospheric condition was monitored with a 19 GHz 
radio seeing monitor using a reference signal from geo- 
stationary satellite. We checked the phase instability by 
the atmosphere with the seeing monitor during this ob- 
servation, and the r.m.s. of 19 GHz phase was sufficiently 
small. The expected typical radio seeing in the 100 GHz 
band was less than at least 2 arcsec during the whole ob- 
serving time, and less than 1 arcsec in the latter half of 
the observation. 

Figure 1 shows the simultaneous light curves of Sgr A* 
at 90 and 102 GHz on 6 April 2005. The IDV of Sgr A* 
is clearly seen at both frequencies. We found that the 
IDV had a rising phase and broad intensity peak. Flux 
density increased from 1.7 to 2.5 Jy at 90 GHz and from 
1.7 to 2.7 Jy at 102 GHz in 1 hr from -18.5 to ~19.5 hr 



UT. The IDV timescale for a two-fold increase is esti- 
mated to be about 2 hour assuming that the increase has 
a constant gradient. The squares in Figure 1 indicate 
the variation of spectral index between these frequencies 
in —5 min intervals. The spectral index is almost flat 
or slightly inverted, although the error in the first half 
is moderately large. Although the inverted spectrum of 
a — 0.3 (S v ex v a ) has been observed in the range from 
several GHz to at least 100 GHz in the quiescent phase 
(e.g., Tsuboi et al. 1999; Melia & Falcke 2001), the spec- 
trum around the intensity peak at 19.5 hr became steep 
inverted up to a — 0.8. A similar inverted spectrum has 
also been observed in the 140 GHz band by using NMA 
(Miyazaki et al. 2004). 

We performed periodicity analyses by using the Lomb- 
Scargle (L-S) method (Press et al. 2007) to search for any 
periodic behavior in the light curve. Figure 2 shows the 
power spectrum of Sgr A* on 6 April 2005. The thick 
curve is the average of power-spectral densities (PSDs) 
at 90 and 102 GHz. There is an excess of PSD around 
the lower frequency limit. The excess is attributed to the 
rising IDV shown in Figure 1. The statistical significance 
of peaks in the L-S periodogram is estimated in terms of 
the false alarm probability (FAP), which is the probability 
for the PSD to exceed a certain level by chance. FAP is 
given as, FAP(> PSD) = 1 - (1 - cxp(-PSD)) M (see 

cq. 13.8.7 in Press et al. (2007)), where M 6.362 + 

1.193 x N + 0.00098 x N 2 for N data points (see eq. 13 in 
Home & Baliunas (1986)), and N is 164. No PSD peak 
with statistical significance was found in the periodicity 
of 1 hour or less. The PSD in the range from / =0.5 to 
2 may show a power law behavior. The straight line in 
Figure 2 indicates the slope with power-law index of / _1 
for comparison. The slope index of the PSD is about — 1. 
Although red noise with a slope index of —2 has been 
reported for Sgr A* in the near-infrared band (e.g., Do et 
al. (2009)), the slope index reported in the radio regime 
is -1 (Mauerhan et al. 2005). 

3. Search for Time Lag 

Cross-correlation function (CCF) analysis is a useful 
technique for finding time lags between light curves at 
different frequencies. However, since it assumes uniform 
sampling, it is not applicable to many astronomical ob- 
servational data. The z-transformed discrete correlation 
function (ZDCF) is a better solution to the problem of 
investigating correlation in unevenly sampled light curves 
(Alexander 1997) and is one method commonly used to 
investigate the time correlation of active galactic nuclei. 
Thus, we used ZDCF to search for a time lag between 
the light curves of Sgr A* at 90 and 102 GHz. Figure 3 
shows the plots of time lag versus ZDCF between 90 and 
102 GHz. ZDCF appears to peak at a time lag of approx- 
imate zero. We estimated the time lag at the peak by 
fitting the data to a quadratic function in the range |time 
lag|< 1 hr. The time lag was —2.56 ±0.92 min. Each 
plot had a scattering within about error around a peak 
of ZDCF. Time lag with actual maximum ZDCF was off- 
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set from a peak of quadratic function shape. We fitted 
the ZDCF with a quadratic function to determine a delay 
against spurious peaks. We estimated also the centroid of 
ZDCF based on all points with correlation coefficients in 
excess of 0.8 r max (r c = ^r.ZDCF^ZDCF,), and the 
centroid time lag r c was —2.57 min. This was consistent 
with the value estimated by quadratic function fitting. 
A positive sign indicates the flux density at 90 GHz has 
a time delay with respect to the one at 102 GHz. The 
observed value indicates no time lag or the flux density 
at 102 GHz being marginally delayed relative to that at 
90 GHz. 

We next examined whether the ZDCF algorithm would 
detect a small time lag in data with the same sampling 
sequence. First, we averaged the light curves at 90 and 
102 GHz to prepare mock-up data, and applied an artifi- 
cial time delay of 10 min to prepare another set of mock- 
up data. We calculated ZDCF between these two data. 
Figure 4a shows the cross correlation function for ZDCF. 
We estimated the time lag at the peak by fitting the mock- 
up data to a quadratic function in the range |time lag - 
10 min|< 1 hr, and the estimated time lag was found to 
be 9.86 ± 0.88 min. The artificial delay of 10 min was 
certainly detected with accuracy of ~ ±0.9 min. To fur- 
ther verify the results, we also calculated ZDCF between 
a 90 GHz light curve that was artificially delayed 10 min 
and the 102 GHz light curve (Figure 4b). By fitting to a 
quadratic function, the time lag at the peak ZDCF was 
estimated as 7.56 ± 0.85 min. As mentioned above, the 
90 GHz light curve was originally found to have a time 
lag of -2.56 ±0.92 min with respect to the 102 GHz light 
curve. Thus, the delay for the mock-up light curves corre- 
sponds to 10.13 ± 1.77 min (= 7.56 min + 2.56 min). The 
time lag between these data is similar to the expected time 
lag. 

In addition, we checked the accuracy of the time lag 
estimated by the ZDCF method by using Monte-Carlo 
simulation. First, we prepared two mock-up data sets by 
the following procedure. An averaged light curve was con- 
structed from the light curves at 90 and 102 GHz. Then, 
the best-fit quadratic polynomial curve was determined 
for the averaged light curve. The residuals of the light 
curves were calculated by subtracting the fitted curve from 
each observed light curve. Two sets of random noise with 
the same standard deviation as the residuals were gener- 
ated. The sets of noise were added on the fitted curve to 
give the mock-up data. Second, we calculated ZDCF be- 
tween these data sets, and we estimated time lag by fitting 
to a quadratic function. Figure 5 shows the histogram of 
the time lag when this process was repeated 1000 times. 
The histogram has an approximately normal distribution 
with an average value of —0.14 min and a standard de- 
viation of 2.69 min. Therefore, if there was a delay of a 
few minutes in the data, we would be able to detect it 
from the observed light curves. Thus, the results shown 
in Figure 3 indicate that there is probably no time lag 
between the light curves of Sgr A* at 90 and 102 GHz. 



4. Discussion 

On the basis of VLA observations of Sgr A*, Yusef- 
Zadch ct al. (2006) reported a time lag of ^20-40 min 
between the flux variability at 22 and 43 GHz. For the 
origin of flare activity, two possibilities have been pro- 
posed: hot plasma in the accretion disk or the expanding 
plasma ejected near GCBH. The time lag that has been 
observed supports the expanding plasma model. Here, fol- 
lowing the analysis by Yusef-Zadeh et al. (2006), we apply 
the expanding plasma model of van der Laan (1966). In 
this model, the synchrotron-emitting plasma blob is opti- 
cally thick when it is ejected, and becomes optically thin 
through adiabatic expansion. When the blob is optically 
thick, the flux density increases as the surface of the blob 
expands. However, the flux density decreases as the opti- 
cal thickness decreases after the transition to being opti- 
cally thin. The frequency at which the light curve is just 
peaking is, 



^o 



R 



1/A 



CD 



where A ~ — (p + 4)/(4p + 6) and p is the index of the 
rclativistic electron energy spectrum [n(E) ocE~ p ] (Yusef- 
Zadeh ct al. 2006). We assume a linear expansion model 
with constant expansion speed, v exp . The radius of the 
expanding plasma is 

R-Ro = v exp (t-t ), (2) 

where t — to is the duration of the expansion. Then, the 
relation between the peak frequency and the duration is 
given by 

A 

(3) 



t-tn = 



-1 



Rq 



J exp 



Therefore, the expected time lag between 102 and 90 GHz 
is calculated as 



Atio2-9o (min) 



102 GHz A -90 GHz' 



43 GHz A - 22 GHz A 

xAT 43 - 2 2 (min), (4) 

where AT43_ 2 2 is the time lag between 43 and 22 GHz in 
minutes. If the time lag /S.T^.22 is 25 min (Yusef-Zadeh 
et al. 2006; Yusef-Zadeh et al. 2008), the time lag Aii 02 _9o 
is expected to be 2.86 min at p = 2 and 2.99 min at p = 3. 
Thus, the flux at 90 GHz is expected to be observed about 
3 min behind that at 102 GHz if the time lag between 
22 and 43 GHz is the same in this observation epoch. 
However, we found no such delay in our data as shown in 
the previous section. The difference between the observed 
time lag of —2.56 ±0.92 min and the expected time lag 
for p = 2 and AT43_22 = 25 min is Atio2-9o.e// = 5.42 min 
(=2.56 min + 2.86 min). If the electron energy spec- 
trum is hard or steep, the expected time lag does not 
change significantly (A£io2-go,e// = 5.19 min at p = 1 and 
Aiio2-90,e// = 5.70 min at p = 5). Thus, the time lag pre- 
dicted by the expanding plasma model is not detected in 
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this observation epoch. There are at least two possible ex- 
planations of this discrepancy. The first possibility is that 
the IDV of Sgr A* is widely diverse, especially in regard 
to the optical thickness of the blob. If the blob is initially 
optically thin even at 100 GHz, the time lag between 90 
and 102 GHz is not necessary. The second possibility is 
that the major portion of the flux above 100 GHz does 
not originate from expanding plasma and instead comes 
from orbiting hot plasma spots on the accretion disk 
around GCBH. The variability may arise from a change 
in the emitting area of optically thick orbiting spots (e.g., 
Broderick & Loeb 2005; Li et al. 2009; Zamaninasab et 
al. 2010). Such emitter would have no time lag depending 
on the frequency. The apparent size of Sgr A* is increas- 
ing according as the well-known law: AD oc A 2 (e.g., Lo 
et al. 1985). If this phenomenon of changing the size is 
caused by electron scattering of accreting matter to the 
GCBH, the size indicate the diameter of the photosphere 
at the frequency. In the case, the lower frequency pho- 
tons from the emitter pass for longer detour in the pho- 
tosphere. The scenario probably explains the higher fre- 
quency photons pass it quickly and have a shorter time 
lag. As other model, the jet-model is proposed by Falcke 
et al. (2009). On the jet- model the time lag depends on 
relativistic speeds of the jet, and the model may explain 
the different time lags. On the other hand, the fine-scale 
flux variations of the observed light curves are very com- 
plicated and can be made of consecutive flares. If so, even 
if the adiabatic expansion is at work for the individual 
flares, the light curves might show no definite time lag 
by blending of the neighboring flares. Regardless, further 
observations with wide frequency coverage are required to 
resolve this issue. 

5. Summary 

Using the NMA, we observed the flux densities of Sgr A* 
at 90 and 102 GHz on 6 April 2005 in order to detect the 
time lag between these frequencies, and we constructed 
light curves covering about 3.5 hour with 1 min bins at 
both frequencies. We estimated ZDCF between the light 
curves at 90 and 102 GHz. The time lag derived from 
the peak ZDCF was —2.56 ±0.92 min. Under the expand- 
ing plasma model, the estimated time lag is sufficiently 
smaller than the time lag expected from the previously 
reported observations at 22 and 43 GHz. In our observa- 
tion data, we did not find a significant delay of the light 
curve at 90 GHz with respect to the one at 102 GHz. This 
result suggests that the plasma blobs from which the IDV 
originates may be diverse, for example, the blob may be 
initially optically thin even at 100 GHz. Another possi- 
bility is that the major portion of the flux above 100 GHz 
does not originate from expanding plasma, and may in- 
stead originates from hot spots on the accretion disk. 
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Fig. 1. Light curves of Sgr A* observed by using NMA at 90 
and 102 GHz on 6 April 2005. Blue and red circles indicate 
flux densities at 90 and 102 GHz, respectively. The integra- 
tion time of each data point is 1 min. IDV, which shows a 
significant increase from UT=18.2 hr to UT=19.8 hr, was de- 
tected at both frequencies. Open squares show the spectral 
index between these frequencies (averaged over 5 min). 
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Fig. 2. Power spectrum density (PSD) of Sgr A* on 6 April 
2005 as computed by the Lomb-Scargle periodgram method. 
The blue and red thin curves correspond to 90 and 102 GHz, 
respectively. The thick curve is the PSD averaged for 90 
and 102 GHz. The straight line indicates slope of / noise. 
Horizontal dashed lines indicate the PSDs that correspond to 
FAPs of 0.01 and 0.001. 
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Fig. 3. Cross-correlation function calculated with ZDCF of 
Sgr A* between 90 and 102 GHz on 6 April 2005. There is 
a correlation function peak at time lag of approximate zero. 
The peak time lag estimated by fitting the data to a quadratic 
function (thick curve) is —2.56 ±0.92 min. Vertical dashed 
line indicates the peak time lag. 




Fig. 4. Cross-correlation functions calculated with ZDCF of 
the mock-up light curve data, (a) Cross-correlation function 
between the average of the observed light curves and the data 
that was artificially delayed 10 min. The time lag, which 
was estimated from the peak of the best-fit quadratic func- 
tion (thick curve), is 9.86 ± 0.88 min. (b) Cross-correlation 
function between the 90 GHz light curve that was artificially 
delayed 10 min and the 102 GHz light curve. The time lag 
estimated from the best-fit quadratic function (thick curve) is 
7.56 ±0.85 min. Vertical dashed line in each panel indicates 
the peak time lag. 
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Fig. 5. Histogram of the ZDCF peak (time lag) for mock- 
-up light curve data with added artificial random noise. The 
number of trials is 1000. The time lag follows a normal dis- 
tribution with an average value of —0.14 ±0.06 min and a 
standard deviation of 2.69 ±0.07 min. 



